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A biometric system of identification and authentication provides automatic 
recognition of an individual based on certain unique features or 
characteristics he or she possesses. Iris recognition is a biometric 
identification method that applies pattern recognition to images of the iris. 
Owing to the unique epigenetic patterns of the iris, iris recognition is 
considered one of the most accurate methods in the field of biometric 
identification. The segmentation algorithm proposed in this article starts with 
determining the regions of the eye using unsupervised neural approach, after 
the outline of the eye is found using the Canny edge, The Hough Transform 
is employed to determine the center and radius of the pupil and the iris. Then 
the normalization allows transforming the segmented circular iris region into 
a fixed-size rectangular shape using Daugman’s rubber sheet model. A 
discrete wavelet transformation (DWT) is applied to the normalized iris to 
lower the size of iris models and to improve classifier accuracy. Finally, the 
URIBIS iris database is used for individual user verification by using the 
KNN classifier or support vector machine (SVM) which based on the 
analysis of iris code as feature extraction is discussed. 
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1. INTRODUCTION 


Nowadays, biometric recognition has become a reliable way to identify and recognize individuals. 
Conventional authentication methods like passwords and identification cards are not always credible because 
these methods can be forgotten or stolen. Identification based on physiological or behavioral characteristics 
of people (such as irises, fingerprints, faces, hand geometry, signatures, fingerprints, et cetera) can accurately 


authenticate people’s identity [1]. 


Iris recognition offers the highest accuracy in identifying individuals compared to any other 
biometric approaches. The iris is defined as the circular part between the pupil and sclera of the eye. 
It is so unique that no two irises are alike, even among identical twins or between the left and right eye of the 
same person [2-3]. These visible characteristics are thought to be discriminatory and unique to each eye, 
as well as stable over an individual’s lifetime. This makes the iris a very useful biometric identifier when it is 
possible to capture iris images of sufficient quality [4]. 


Journal homepage: http://ijai.iaescore.com 


430 im ISSN: 2252-8938 


The system is inspired by Daugman’s work [5]. It is composed of a number of sub-systems, 
which are described in Figure 1 and Figure 2 the contains stages of an iris detection system. 

The segmentation step enables the localization of the iris region in an eye image, which consists of 
finding two contours of the iris (pupil and iris). In other words, the internal contour (pupil/iris) and outer 
(iris/sclera) are located. Once the iris region is segmented, the next stage is to normalize this segmentation. 
The normalization transforms the iris texture from Cartesian to polar coordinates in order to limit the 
problems of variations in the eye (like the optical size of the iris or the position of the pupil in the iris). 
This step is based on Daugman's rubber sheet model [5]. Feature extraction is a technique used to extract 
information from the iris image, in order to lower the size of iris models and improve classifier accuracy. 
These features can not be used to reconstruct images, but these values are used in classification. 
Finally, matching is performed by comparing the features of a template iris with the feature vectors of 
templates in the database, and a decision is formulated [6-7]. 
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Figure 1. Iris recognition system Figure 2. Typical stages of iris recognition 


The main motivation in this research is to propose a new effective and robust algorithm 
to segment and classification the clear or noisy iris images. The proposed algorithm added a new 
pre-processing step using a new unsupervised neural approach to divide the iris image in two regions namely 
iris and eyelashes region, sclera and skin, in order to facilitate the determination of the iris contour in the 
phase Canny edge detection [8]. Unlike other iris segmentation algorithms, that process the whole image of 
the eye (which contains the non-iris regions that generate segmentation errors) such as Daugman algorithms. 
Then, discrete wavelet transformation (DWT) is used for extracting the optimum features of iris images and 
reducing the runtime of classification of these iris templates. Finally, a comparison was made between the 
classifiers SVM and KNN. 


2. METHODOLOGY 
The iris recognition system consists of the following steps: 
— Image acquisition 
— Iris segmentation 
— Normalization 
— Feature extraction 


— Classification 
In the approach presented here, segmentation was achieved using the Hough transform with the aim 
of locating the iris region of an eye image. The normalization enabled the transformation of the segmented 
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circular iris region into a fixed-size rectangular shape using Daugman’s rubber sheet model. 
Two level discrete wavelet transformation (DWT) was applied to the normalized iris for feature extraction. 
The support vector machine was used to classify the similarity between the iris templates. Figure 3 shows a 
diagram of the iris recognition system. 
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Figure 3. Diagram of iris recognition system 


2.1. Iris segmentation 
In this section, we present a new iris segmentation approach used as a pre-processing step, based on 


neural competitive concepts [9]. As shown in Figure 4, this approach allows partitioning the iris image in two 
regions namely iris and eyelashes region, sclera and skin. Then, the outline of the eye is found using the 
Canny edge, The Hough Transform is employed to determine the center and radius of the pupil and the iris. 
Figure 4 shows the iris segmentation steps. 
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Figure 4. Algorithm used for iris recognition 
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2.1.1. Iris segmentation using the unsupervised neural approach 

In an unsupervised context, that means when no prior information about the data sample, 
we involve the neural competitive clustering procedure to the Iris segmentation. The objective of these 
approaches to classifying the image pixels according to their distribution in the representation space and to 
assign them a label. As illustrated in Figure 5, it starts first by the organizationIt starts first, by the 
organization of the pixels of Iris image in an observation matrix (each row represents a pixel and each 
column represents an attribute) to estimate the underlying probability density function (pdf)of the pixels 
distribution using a non-parametric estimator. In the second step, the procedure uses an artificial neural 
network with competitive training (NNCT) to extract the local maxima of the pdf. Following a modes 
detection method using a technique to detect the existing interneural connection [10]. The last step is for 
affecting the remaining pixels to their classes. 
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Figure 5. Architecture of the segmentation procedure 


2.1.2. The estimation of underlying probability density function 

After constructing the observation matrix of an Iris image pixels’ = {Xi Xo; “es _Xq considering as a 
set of Q N-dimensional observations (in this case N=3) With Xq = Eee Xq,20 +» Xqny aston q=1,2,....Q 
and a probability density function P(X), anon-parametric method based on The Parzen [11] window using to 
estimate this underlying density function. The technique is a fast estimation algorithm that is proposed by 
Postaire and Vasseur [12]. First, the range of variation of each component of these observations is normalized 
to the interval [0, R], where R is an integer such as R = 2, by means of the transformation defined as: 


(2n,q-min tna) 


Ynq = *R () 
(max xng—minxn,q 

Each axis of the so normalized data space is then partitioned into R exclusive intervals of unit width. 
This discretization defines a set of RNhypercube of unit side length. Each hypercube noted H(X), 
is a site defined by its N coordinates x1,%2,...,Xy,)--,X\y Which are the integer parts of the coordinates 
of its centre X. To be more specific, let Yy = [y1.q+¥2,q> Vago Ynql> q=1, 2, Q be the Q observations in 
the normalized space. Each observation Yq is found inside a non-empty hypercube with the coordinates 
Xn = int (Ynq): n=1, 2, N, where int(Ynq) designates the integer parts of yp,q. If several observations fall in 
the same hypercube, this one appears many times on the list of non-empty hypercubes. Furthermore, 
the number of times the hypercube H(X) appears in that list indicates the number of data points q[/H(X)] 
which falls in this hypercube. Subsequently, the value of the local density estimated is: 


p(X) = we (2) 


Since the volume of H(X) is equal to unity. 

So, this fast procedure allows only the estimation of the underlying probability density function at the centers 
of the non-empty hypercube whose number never exceeds the number Q of available observations (pixels). 
At the centers of the hypercube cells, which are not on that list, the density estimates are known to be null. 
At the end of this fast algorithm, all the available information for clustering is in the discrete set X of 
estimated values of the underlying probability density function P(X). 
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2.1.3. The extraction of local maxima by neural network 

Assimilating the modes to the local maxima of the pdf, the proposed approach uses 
the Neural Networks with Competitive Training (NNCT) [13]. In the training algorithm, we work only on the 
pdf by presenting sequentially, the centers of the non-empty hypercube of the set X to the network, instead of 
the Q observations. 

The neural network is composed of two layers: the input layer and the output layer. The first one is 
made of N units [,,, n=1,2,...,.N, such that unit J, is solicited by the attribute X, of the non-empty hypercube 
H(X) when this one is presented to the network. However, each output neuron materializes a hypercube 
which represents the site of one local maximum of the pdf in the set X, and presents its weight by the mean 
vector U;(X), k=1,2,...,K. The number of the output neuron is first initialized arbitrarily. 

During the training phase, the output neurons enter into competition with each other by comparing 
the distance D[u;,(X) , H(X)], k=1,2,...,K, between the input hypercube H(X) and each output neuron pu, (X), 
the winner is the closest one to the hypercube, then we compare the values of the pdf associated to the winner 
neuron [Ug (X) and to H(X). The distance measure used in this training algorithm is Mahalanobis distance that 
gives the best results for the non-Gaussian distribution [14] instead of Euclidian distance as in the (NNCT). 


2.1.4. Detection of significant modes of pdf 

The aim of this phase is to connect each group of the closest modes in such a way that we get a map 
which preserves the shape and structure of the classes existing in the image, by an improved Competitive 
Hebbian Learning method (CHLim) [10] which allows eliminating the influence of the output neural number. 

To generate the induced Delaunay triangulation, the CHLim, given the K modes detected by CNN 
as prototypes in RN, successively adds connections among them. The method does not change the weight of 
prototypes, but only generates topology according to these prototypes. For each mode H(X;,), 
can be connected to the two closest prototypes by an edge using Mahalanobis distance as a measure of 
resemblance instead of Euclidian distance, it works as an activation function for competition 
between neurons. This leads to the induced Delaunay triangulation, which is limited to those regions of the 
input space R™. 


2.1.5. Classification procedure 

Once the different modes are identified, the classification method [15] that we use in the present 
work consists first, in defining the pixels falling into any mode of a connected set as the prototype of one 
cluster. Then, the remaining pixels, which do not fall in one of the detected modes, are assigned to the 
clusters attached to their nearest neighbor among these prototypes by means of Mahalanobis distance. 
Figure 6 presents the procedure for classifying the iris image. 
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Figure 6. Illustration of the Iris segmentation procedure 
2.1.6. Hough transform 


Hough transform is a standard image analysis tool for finding curves that can be defined in a 
parametrical form such as lines, circles, parabolas, and hyperbolas [16]. The edge map is then used in a 
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voting process to maximize the defined Hough transform for the desired contour. Considering the obtained 


edge points as (xj ; yj), J =1523...... , n, a Hough transform can be written as: 
H (Xe, Ver) = pas R(X}, Vys Xoo Ve sT) (3) 
Where 


1 i Xi VirXer r)= 0 
Stayer) = { f g( ‘j Vj Cc Ye ) _ 
0 otherwise 
The limbus and pupil are both modeled as circles and the parametric function g is defined as: 


I (Xj Vj XerVor?) = (Gj — Xe)? + OY — Xe)? = 7? (5) 


Assuming a circle with center (xc; yc) and radius r, the edge points that are located over the circle 
result in a zero value for the function. The value of g is then transformed to 1 by the h function, 
which represents the local pattern of the contour. The local patterns are then used in a voting procedure 
using the Hough transform, H, in order to locate the proper pupil and limbus boundaries. In order to detect 
the limbus, only vertical edge information is used. The two eyelids usually cover the upper and lower parts; 
these parts contain horizontal edge information. The horizontal edge information is used for detecting the 
upper and lower eyelids, which are modeled as parabolic arcs [9]. 

After the Hough Transform process is complete, six parameters are stored: the radius and x and y 
(the center coordinates for both circles). Eyelids are isolated by first fitting a line to the upper and lower 
eyelid using the linear Hough transform. A second horizontal line is then drawn that intersects with the first 
line at the iris edge that is closest to the pupil. This process is completed for both the top and bottom eyelids. 
The second horizontal line allows maximum isolation of eyelid regions. Figure 7 shows an example of 
properly segmented iris. 


Figure 7. Examples of the segmented iris 


2.2. Iris normalization 

The size of the iris may change because of variation in the illumination, even for an iris from the 
same person. This elastic deformation in iris texture must be compensated for to achieve more accurate 
recognition results. Normalization is a technique to prepare a segmented iris for feature extraction and 
transforms the segmented circular iris region into a fixed-size rectangular shape. The Cartesian-to-polar 
transform of the iris region is based on Daugman’s rubber sheet model [5], as illustrated in Figure 8. 
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Figure 8. Cartesian-to-polar coordinate transform 


The Cartesian-to-polar transform of the iris region is modeled by (4). 


I(x(r,@)), v(r,8)) > I(r, 8) (4) 
with 


x(r,6) = (1 —1r)xp(@) + rx;() 
y(r, 6) = (1 -1r)y,(@) + ryi(8) (5) 


The coordinates of the inner and outer boundaries along the @ direction are xp, yp and 
Xi, yi, respectively. Image-normalized enhancement is required to obtain accurate features because the 
normalized iris presents many variations due to variations in light illuminations. By using histogram 
equalization, pixels intensities are spread over the entire iris template. Figure 9 shows an iris template with 
histogram equalization. 
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Figure 9. (a) Iris normalized into polar coordinates, (b) Enhancement of the iris normalized image 


2.3. Feature extraction 

In the encoding stage, two level discrete wavelet transformation (DWT) is applied to the normalized 
iris region, in order to determine DWT coefficients, by passing the signal in different frequency ranges, 
namely low-low (LL), low-high (LH), high-low (HL) and high-high (HH), as indicated in Figure 10. 
The frequency range can be represented as LL < LH < HL < HH [17-18]. 

The LL sub-band represents the feature characteristics of the iris; this sub-band (LL) is used in the 
classification phase. Figure 11 indicates the resolution of the normalized and enhanced iris image (48x432). 
After applying DWT to an enhanced iris image, the resolution of sub-band is (24 x 216). This sub-band is 
used instead of the original normalized iris data to reduce the resolution of the iris model, and the runtime of 
the classification is similarly reduced. 
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Figure 11. (a) Enhancement of the iris 
normalized image, (b) DWT transforms normalized 
iris image into LL sub-bands 


Figure 10. Iris template using DWT 


2.4. Support vector machine SVM 

The support vector machine (SVM) performs pattern recognition according to the principle of 
structural risk minimization. The development of SVM as a classifier has two major aspects: the first aspect 
is to find the hyperplane that optimally separates the two classes, and the second aspect is the transformation 
of the linearly non-separable classification problem into a linearly separable problem [19]. 

Given an input entity vector learning set and the class tag (x_i,y_ i), i=1,2...,1, where x,€ R” 
and y € {1, -1}. In the case of a linearly separable problem, the separating hyperplane (which defines the 
boundary between class | and class 2) can be represented as follows: 


MIN yb, To tee is (6) 
Subject to 
yi(w"' p(x) +b) 21-§ 55 = 0 (7) 


In the case of a non-linearly separable problem, the training vectors x_i are mapped into higher dimensional 
space by the function @ in order to make the problem linearly separable. The SVM then finds a linear 
separating hyperplane with the maximal margin in this higher dimensional space. The penalty parameter of 
the error term is C > 0. Furthermore, K (xi, x;) = @ (x;)" (x;) is called the kernel function [20]. 
— Linear: 

K(x;,x;) = x7 xj (8) 
— Polynomial: 


K(x, x)= (yxi7 45 + 1)7,y > 0 (9) 


- Gaussian: 


K(x;,x;) = exp(- ey y>0 (10) 


3. RESULTS 

The Classification Learner App (included in the Statistics and Machine Learning Toolbox 
for MATLAB) is an application that can be used to train models to classify data using supervised 
machine learning methods. The SVM and KNN classifiers were used for our approach [20-21]. 
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By using Classification Learner, we can explore data, select features, specify validation schemes, 
and train and export classification models to the MATLAB workspace to generate predictions for new data. 
To evaluate the effectiveness of the proposed system, a set of eye images obtained from the UBIRIS 
databases were used [22]. Figure 12 shows an example of the eye images contained in this database. The total 
number of iris classes is 100, and each iris class contains three images stored in the database for this 
experiment. In this work, 300 images were used for training, and 100 images were used for testing purposes. 


Figure 12. Eye images in the UBIRIS database 


Four kernel functions (Linear, Quadratic, Cubic, and Medium Gaussian) were used, and the 
classification efficiency of each kernel was measured. Since the highest classification accuracy was obtained 
by a Medium Gaussian kernel, as shown in Table 1, this kernel was used in our system for classification and 
recognition purposes. As the simulation result shown in Table 2, the KNN algorithm assured better the 
accuracy of 95%, as well as the execution time is too low compared to that of the SVM algorithm. 


Table 1. Efficiency of various kernel Table 2. Efficiency of KNN classifier 
functions (SVM) Classification Accuracy Training Time 

Kernel Type Classification Training Time (%) (sec) 

Accuracy (%) (sec) KNN 95.00 19.94 
Linear 68.2 430.89 
Quadratic 84.3 533.09 
Cubic 88.6 580.83 
Medium Gaussian 88.9 602.28 


The accuracy of our proposed technique is 95%, is shown in Table 3, which is better than the other methods 
that are used for runtime comparison. 


Table 3. Accuracy comparison of our proposed technique with others 


Author Classification Accuracy (%) 
H Ohmaid et al. (Our proposed) 95 
MH Hamd and al. [23] 94 
SG Firake and al. [24] 92.85 
J poonia and al. [25] 82.5 


4. CONCLUSIONS 

In this paper, segmentation was achieved using pretreatment step based on an unsupervised neural 
approach. The normalization enabled the transformation of the segmented circular iris region into a fixed-size 
rectangular shape using Daugman’s rubber sheet model. Discrete wavelet transformation (DWT) is used for 
extracting the optimum features of iris images and reducing the runtime of classification of these iris 
templates. Finally, a comparison was made between the classifiers SVM and KNN. Experimental evaluation 
using the URIBIS database clearly demonstrated the proposed technique performs in an efficient manner, 
especially for the KNN classifier that gave an accuracy of 95%. 
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